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Abstract 

A simple and general approach for calculating the elliptic finite-gap solutions of the 
Korteweg-de Vries (KdV) equation is proposed. Our approach is based on the use of 
the finite-gap equations and the general representation of these solutions in the form 
of rational functions of the elliptic Weierstrass function. The calculation of initial 
elliptic finite-gap solutions is reduced to the solution of the finite-band equations with 
respect to the parameters of the representation. The time evolution of these solutions 
is described via the dynamic equations of their poles, integrated with the help of the 
finite-gap equations. The proposed approach is applied by calculating the elliptic 1-, 
2- and 3-gap solutions of the KdV equations. 

1 Introduction 

In accordance with the finite-band theory, the integrable Korteweg-de Vries (KdV) equa- 
tion can be considered as the compatibility condition of the two auxiliary linear differential 
matrix equations 9x±^ = U± (j;+ = x, x_ = t, dx = d/dx) with matrix operators U±, 
coefficients of which, as is well known (see [||, |2|, |3|), are expressed through solutions of 
the KdV equation. The finite-gap solutions of the KdV equation are solutions of the 
spectral problem for these equations with a finite-gap spectrum of eigenvalues. These are 
expressed through multi-dimensional Riemann theta functions with implicit parameters, 
the evaluation of which is the special algebraic geometrical problem (see @, ^)- The 
class of elliptic finite-gap solutions leads to the problem of the reduction of n-dimensional 
theta- functions to the one-dimensional Jacobi theta functions (see @). 

In the framework of the spectral problem the finite-gap solutions satisfy the matrix 
finite-gap equation of the form d±'^ = \J± ^ — ^\]± , where ^ is the matrix components 
which are polynomial in the eigenvalue of the auxiliary equations (such as in the case of 
the "sine-Gordon" equation Q). Usually, solving this equation is realized with help of the 
well known Abel transformation with a subsequent solving of the inverse Jacobi problem 
(see 1,1). 

However, in the case under consideration the finite-gap solutions of the KdV equation 
are elliptic functions represented as rational functions of the elliptic Weierstrass function 
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(p-function) ||8|. Solving the finite-gap equations in terms of p-function can be reduced to 
solving simple algebraic equations. This yields the straightforward manner for calculating 
the elliptic finite-gap solutions in an explicit form. 

We shown that in the initial time {t = 0) the finite-gap equations can be reduced to 
algebraic equations with respect to parameters of the above mentioned rational functions. 
The computation of these parameters gives all possible initial elliptic finite-gap solutions 
in the form of linear combinations of p- functions with shifted arguments. 

In accordance with the KdV equation the time dependent elliptic solutions are built as 
linear combinations of p-functions with time dependent argument shifts (ipi) (its poles) 
under condition of its compatibility with corresponding initial solutions. Their time evo- 
lution is determined by the poles which satisfy the system of linked dynamic equations 
which follow from the above mentioned auxiliary equations. 

Using the finite-gap equations we shown that the dynamic system is transformed in the 
system of independent differential equations of the first order with separated variables of 
the form dt(pi = Xi{(pi). Here Xj-functions represent themselves roots of some polynomial 
equations which follow from the finite-gap equations the order of which equals the number 
of poles in the elliptic solutions. 

This paper is organized as follows. In Section § the approach to the straightforward 
calculation of the elliptic finite-gap solutions of the KdV equation, based on the auxiliary 
system of finite-gap and dynamic equations, is formulated. In Section ^ this approach is 
applied to the calculation of the initial elliptic 1-, 2- and 3-gap solutions of the KdV equa- 
tion. In Section ^ the time evolution of these elliptic solutions is investigated. We shown 
that the linked system of auxiliary dynamic equations, for poles of the corresponding time 
dependent elliptic solutions, can be integrated with the help of the finite-gap equations. 

2 Finite-band equations and general elliptic solutions 

In the class of elliptic functions (which we shall denote as U{x,t)) with the p- functional 
representation under consideration, the finite-gap solutions of the KdV equation can be 
considered as solutions of the finite-band equations. The latter gives the compatibility 
condition of the finite-gap and the general solutions of the auxiliary linear differential equa- 
tions. These finite-band equations represent a system of equalities, obtained by equating 
coefficients of the power series in the eigenvalue E of the finite-gap and general solutions 



The finite band equations (which are also known as "trace- formulae" ) can be written 
in the form g || 



(g is the number of a gaps in the spectrum of the eigenvalues E), where the Xn-functions 



(see 



§)• 




(2.1) 
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are determined by the recursion relation 

n-l 

Xn+1 = dxXn + X] XkXn-k, Xi = -U {x)- (2.3) 
k=l 

Here An and Xn are coefficient functions of the power-series expansion in E of the general 
and the finite-gap solutions of the auxiliary equations, respectively. In view of ( p. 3D the 
X-functions have the form of polynomials in the solutions U{x,t) and its derivatives. As is 
well known, is algebraicly solvable with respect to the coefficient functions bn{x,t). 
Therefore, the finite-gap equations ( |2.1D are reduced to the system containing expressions 
for the coefficient functions 6„(x, t), n = (1, . . . ,g), which at n < (7 is a closed system of 
equations for the elliptic finite- gap [/-solutions. 

The elliptic finite-gap solutions as double periodic functions of the complex variable z 
which admit the rational functions representation namely the elliptic Weierstrass function 
p . In view of the equation (0) , U is determined by the formula 

with poles of the secondary order in p = p{z\uj,u') (u and lo' are real and imaginary 
half-periods of the p-function), where R{z,t) means an odd function of z. 

At the initial time {t = 0) the expression ( |2.4D describes all possible so-called initial el- 
liptic finite-gap solutions, which in the case of even functions {R = 0) under consideration, 
takes the form 

t/w = «pM + EE 

that reduces the finite-gap equations ( ^ ) to simple algebraic equations with respect to 
the a- and /ij-parameters. 

The time evolution of the elliptic finite-gap solutions ( |2.4| ) are determined by the time 
dependence of their poles, which are described with the help of the known dynamic aux- 
iliary equation 

dt'^ix, t, E) = (4a + 2,{Ud^ + d.:,U)) ^'(x, t, E), (2.6) 

^ = ^'({a„)}, {bn}), bn = bn{U{z,t), U^"'\z,t)) can be reduced to the system of dynamic 
equations with respect to the poles of the [/-functions. This is a system of coupled dif- 
ferential equations of the first order which, as will be shown below, can be reduced to 
independent equations with separated variables. 

3 A calculation of the initial elliptic solutions 

The calculation of the initial finite-gap elliptic solutions of the KdV equation is based on 
the use of the finite-band equations ( |2.lD in the representation of the rational functions 
(|2.5| ). On equating the corresponding coefficients of the Laurent expansion in p of the 
left- and right-hand sides ( |2.1| ), yields simple algebraic equations which determine the 
parameters of the expression ( |2.5D for the initial elliptic finite-gap solutions of the KdV 
equation. 
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3.1 One-gap elliptic solutions 

In the one-gap case the parameters of the elhptic solution U{z) are determined by the 
system of three finite-gap equations of the form ( |2.1| ) at n = 0, 2. A substitution in these 
equations, of the explicit expressions (|2.2|) for An and polynomials in U expressions, for 
the Xn-functions which follow from yields the system 

^{(«2 - U) + 2{bl - b2) - ai6i} = -l{C/2 _ 

2 4 2-^ /o i\ 

1 3 3 3 1 v*-"-^; 

^{(gffl? - -aia2 + -as) + 3{-aj - 02)61 + 3ai(6f - 62) 

+ (126162 - 4!63 - 66?)} = - ^U^^^' + QUU^^^ - 2U^} 

2^ 

in which 6„j„>2 = (in view of the relation bn\n>g+i = 0), where g is the number of gaps 
in the spectrum of the auxiliary linear differential equation. Excluding 6„ from the system 
(3.1) we can obtain the equations 

62 = = ^(3[/2 - C/(2)) + ^aiU + - ^aj; 

63 = = - -^(C/(4) + lOC/3 - 5C/(i)' - 10[/[/(2)) - -^ai(3[/2 - C/(2)) (3.2) 

1 rr. 2 . x 1 1 1 3 

+ Ye ~ 2"^ 4"^"^ ~ 16^^' 

which is known Q as the one-gap "trace formulae" . These equations form a closed system 
which determines the initial elliptic solutions (2.5). Inserting the rational expression for 



U{z) (2^) into the system ( |3.2D and equating coefficients of the Laurent expansion in p 
in the right-hand sides to zero, we obtain algebraic relations which lead to the equalities 

1) a = 2ai^ = a2, = 0, ai = 0, 02 = -^52^, as = -jds; 

2) a = 2, a(i,2). = ^/3(i,2),ai = 0, a2 = ^(ll52-120p2((^i)), (3.3) 
«3 = -^93 + 4p((pi)(12p2((^i) _ ^2) - 4p'2((^i). 

Here and below (3i = 12h? —92,^2 = 2/i'^, cpi is the argument of the function h = p{^i) 
which satisfies the equation (12/i^ —92)^ = 48/i/i'^, from which tpi = (2/3)a'j|j=(i^2,4)) where 
ooi = [uj — uj') (see also pO|]). 

The two systems (|3.3[) determine two types of initial elliptic solutions: 



1) Uiz) = 2pizy, 

2) U{z) = 2p{z) + 2[p{z - (^i) + p{z + (^i)] - 4p(<^i), 

where the first type is the known Lame potential and the second type is a new potential 
obtained in [IToll. 
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3.2 2-gap elliptic initial solutions 

Coefficients of the initial elliptic 2-gap solutions are determined by the system of the 
four finite-gap equations of the form ( |2.1D at n = 0, 3. In analogy to the one-gap case, 
the explicit form can be obtained by substituting the expressions for and the 

expressions for Xn (from ( |2.3D ) into ( |2.1| ). In doing so, the first two equations, which 
coincide with the first two equations of the system (|3.2| ), are solvable with respect to bi 
and 62- Excluding the latter from the fourth and fifth equation and taking into account 
the equality bn\n>3 = we obtain the finite- gap system 

63 = = ^(16a3 + Saaf/ + lOC/^ - Sf/'^ - 2aiU" 

- WUU" + U^^y); 

64 = = ^(-16a^ + 6A*a4 + 32a3U + 2402^7^ + 35U^ (^-^^ 

- lOUU'"^ - Sa2U" - 70U^U" + 21U"^ 

Inserting the rational expression ( |2.5| ) for U{z) into the system (3.5) and equating coeffi- 
cients of the Laurent expansion in p, on the right- and left-hand sides, we obtain algebraic 
relations which lead to the equalities 

1) a = 6, ai, = Q2, = 0, ai = 0, 02 = -"^§2, 

27 27 2 81 

03 = -t53, «4 = —92, 04 = -Xf253; 

2) a = 6, ai, = 5i^j{l2e] - igs), 02 = -7(9e2 + 2Aj), 

03 = 18(3e] - 5ejAj), 04 = 27(36e| + 16e|Aj + 3A^), 
as = -54(36e| - 52e|Aj - QejA^), 

kj = Se] - \g2; 

3) a = 6, ai, = (5„ +5,,fc)(12e2-(72), a2 = 161Aj- 378e|, (3.6) 
as = 531ejAj 108e^, 04 = 27(240A^ + 1280e| - 1159e^Aj), 

as = 27(1594ejA| 8100ef + 120e|A|); 

4) a = 6, ai^ = \f3i5i^i, 02, = 1^2^1,1, 

a2 = |(18/i2 + |ai), as = 1(24/1^ - 25/i'2) - ^ai), 

04 = f (144/i^ + 44/i/i'2 + 56/i2qi + fa?), 

as = -^(144/1^ + Gh'^h'^ - 74h^ai - fh'^ai - ^aj), 



where ip2 is determined by the equality h = p{ip2) which satisfies the equation 64/i'^ -|- 
48hh''^ai—a'l = 0. These four types of equalities ( |3.6| ) lead to the following four expressions 

1) U{z)=Qp{z); 

2) U{z)=Qp{z) + 2p{z + u:i)-2ef, 

3) U{z)=Qp{z) + 2p{z + io.{) + 2p{z + oJk)-2{p{uJi) + p{ojj))- 

4) U {z) = Qp{z) + 2p{z - ^2) + 2p{z + ^2)) - Mf2)- 



Finite-Gap Solutions of the KdV Equation 



27 



which describe all possible initial two-gap solutions of the KdV equations (see |]ll|, |10|). 
The first type coincides with the known two-gap Lame potential, while the second and 
third type corresponds to the Treibich-Verdier potential . The fourth solution is the 

new two-gap potential obtained in In 



3.3 Initial elliptic three-gap solutions 



The parameters of the initial elliptic three-gap solutions (2.5) are determined by the system 



of the finite-band equations ( ^.l]) at bn<3 7^ and 6„,>3 = 0. Using the expressions (2.2) 



and (2.3) this system can be expressed through [/-functions. In doing so, the first three 



equations coinciding formally with ( p.l| ) and are solvable with respect to 63. Therefore, 
taking into account the equality &n|„>4 = 0, we can obtain the equation 

64 = = - 1602 + 64o4 + 32a3U + 24a2U^ + 35*U^ 
- 70 * UU''^ - 8 * a2U" - 70 * U^U" 

+ 21C/"2 + 28U'U"' + 14C/[/(^) - , (3.8) 



which is solvable with respect to parameters of the rational expression U{z) ( |2.5D . Substi- 
tuting (|2.5| ) into the equation (|3.8|) and equating coefficients of the Laurent expansion in p 
with its right-hand side to zero we can obtain closed algebraic relations for a-parameters 
of the [/-solution (^^ . 

It can be shown that there are four possible types of solutions for the a-parameters. 
The corresponding four types of initial elliptic solutions can be written in the form 



(3.9) 



1) [/(z) = 12p(z); 

2) U{z) = I2p{z) + 2p{z + LUi) - 2ei, ei = p{LOi); 

3) U{z) = I2p{z) + 2{p{z + LOi) + p{z + ujj)) - 2{ei + ej); 

4) Uiz) = \2p{z) + 2(p(z + (^3) + p{z - V93)) - 4p(v?3). 

Here the argument (^93 is determined by the equation 

,fi 101 ,4 29 o 43 2,2 23 , / I o 5 2\ 

/i H Q2h H 93/1 92173/1 - 99 H qi = 

196^ 49^ 784^2 196^^ V3136^^ 98^V 

where h = p{(p3). The initial solutions 1 and 2, 3 are the 3-gap Lame (M) and generated 



on the 3-gap case Treibich-Verdier |12] potentials, respectively. The initial elliptic solution 
4 is the 3-gap generalization of the solution obtained in jlO[ |. 

Note that the above described algorithm is general and applicable to computing arbi- 
trary initial elliptic n-gap solutions of the KdV equation. 



4 A dynamics of the elliptic finite-gap solutions 

The time dependent elliptic finite-gap solutions of the KdV equation have the general form 
of the rational functions of the p-function (|2.4| ) , parameters of which are functions of time 
t. These parameters are described by the system of the auxiliary dynamic equation ( |2.6D 
and the finite-gap equation (|2.1|) . 
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By substituting the expression (|2.4D for the ehiptic finite-gap solutions U{z, t) into the 
KdV equation dtU{z,t) = 6U{z,t)U' {z,t) — U"'{z,t) and equating the coefficients of the 
Laurent expansion in p, in its left- and right-hand sides, lead to the known general formula 



N 



U{z,t) = 2^p{z-^^{t)) + C, (4.1) 
1=1 

in which the integer number N and the constant C are determined by the condition of the 
reduction of U{z, t) to the corresponding initial elliptic finite-gap solution at t — > 0. In the 
cases of the elliptic 1-, 2- and 3-gap solutions the numbers must provide the reduction 
of the time dependent solutions ( |4.1| ) to the corresponding initial elliptic solutions of the 
systems (pi), (|3.7[) and (|3.9D, respectively. 



The substitution of the expression (4.1) into the KdV equation reduces the problem of 
the time evolution of U{z, t) to the time evolution its poles. The latter is described by the 
system 



N-l 



dm{t) = -12X,{t) + C, Xi{t) = p{ipi{t) - ^j{t)) {g > 2), 

N 



N ^ ' 



n=l 

which is a system of coupled equations. However, in view of the symmetry properties of 



the finite-gap equations, ( |2.1| ) determine Xi(t) as function Aj((^j(t)). Then the system 

-12t, (4.3) 



can be transformed to 



X,{hi{^i)) + C 



which describe the time dynamics of the poles in the expression ( |4.1| ). Here, initial values 
(fQi = (pi{Q) are determined from the expressions for initial elliptic finite-gap solutions. 
The functions Xi[(p,i) are determined by the finite-gap equations (|2.lD with ?7-functions 



in the form (4.1) as roots of A^th order polynomials in Xj. These polynomials are followed 
from the algebraic equations which can be obtained by equating coefficients of the Laurent 
expansion in p of the left- and right-hand sides of the indicated equations ( |2.1| ). 

Thus, the problem of time evolution of the elliptic finite-gap solutions of the KdV 
equations is reduced to the solution of the finite-band equations, with respect to the Aj- 
functions from relations ([4.3|). The proposed approach will now be applied in calculating 
the time evolution of elliptic 1-, 2- and 3-gap solutions. 



4.1 Elliptic 1-gap solutions 

The types of time dependent elliptic 1-gap solutions of the KdV equation are determined 
by the values of the number N in the expression (4J.). The condition of a coincidence in 
the general expression ( ^.1| ) with the initial elliptic 1-gap solutions ( |3.4| ) at t — > 0, yields 
the values, namely A = 1 and A = 3. The corresponding elliptic 1-gap solutions have the 
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form 



1) Uiz,t)=2piz-ipiit)), 



(4.4) 



2) =2^p(z-^i(t))-4p((^i). 
1 

The time dependence of poles of the elliptic solutions 1 and 2 is determined by the second 



equation of the system (O) at = 1 and N = 2, respectively. In accordance with the 



initial conditions defined by the system (3.4), the solutions of the dynamic equation have 
the form ipi = cit at = 1 and 

ipi{t) = Cit, Lp2{t) = C2t + ipf\ ifsit) = est - ipf^ 



at N = 3. The substitution of these expressions into (|4.4|) yields the following two expres- 
sions 

1) U{z,t) =2p{z-cit); 

2) U{z,t) = 2{p{z -cit) + p{z - C2t + ^) + p{z - est - 99^)) - 4p(^i), 

which describes two possible types of the elliptic 1-gap solutions in the form of superpo- 
sitions of one and three independent traveling waves, respectively. 



4.2 Elliptic 2-gap solutions 

The possible types of the time dependent elliptic 2-gap solutions of the KdV equation ( [4. 1]) 
are determined by the values which can be obtained from the compatibility condition 
between (4.1) and ( |3.7D as t ^ 0. Under this condition the number equal 3, 4 and 5 in 



the formula (4.1). 

1. The time dependent 2-gap elliptic solution corresponding to the initial condition 1 



in the system (|3.7D, which is described by the expression (4.1) at A^ = 3, has the form 



Uiz,t) = 2^piz-^iit)). 



(4.5) 



i=l 



The time evolution of poles (/'i(^) is described by the equation ([4. 3]) in which Xi = 
Y^j^i j=i Pi'^i ~ Vj)- Under initial conditions, the lower limits of the integration in ( ^.3[ 
are ipoi = 0. 

To comput the Xj-function we substitute the expression ( |4.5D into the finite-gap equa- 



tion (3.2). Then, equating coefficients of the Laurent expansion in p of the left- and 
right-hand sides, we obtain algebraic equations which are reduced to the polynomial equa- 
tion 



+ C2X^ + CiX + Co 



0, 



(4.6) 



where Cn = c„((/?), n = 1,3 (here and below the subscript i of 99 is omitted). Three 
solutions of (|4.6|) describe three unknown functions Xi{ipi), i = 1,3. The dependence of 
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the coefficient functions c„ on Lp is described by the expressions 

CO = -^53 - Y|a2 - ^(14952 + na,)h + ^fe + Aa,) (| 
29 8 

42 3 

Co = h H — 

125 800 V h' 

where f3i = 12/i^ — 172 and h = p{(p) so that Xi depend on ip through the /i-function. 

2. The time dependent elhptic 2-gap solution corresponding to the condition 2 in the 
system ( |3.7| ) is described by the expression (4.1) at = 4 and has the form 

4 

U{z, t) = 2 p(z - if.it)) - 2p((^o4). (4.7) 

i=l 

The time evolution of poles ipi{t) are described via the function Xi = Yl'j^ i i=i ~ ^j)- 
In view of the initial conditions, the lower limits of an integration in ( |4.3| ) are (/^oj|j<3 = 
and (/9o4 = ^j- 

For computing the Aj-functions we substitute the expression ( [4. 7] ) into the system (3.5). 
Then, equating coefficients of the Laurent expansion of the left- and right-hand sides, we 
obtain the equation 

X^ + 03^3 + C2X^ + CiX + Co = 0, (4.8) 

solutions of which describe the functions Xn{ip), n = 1,4. The dependence of the coeffi- 
cient functions c„ on (/? are described by the expressions 

Co = — fh\g2 + ?7i25202 + "1302 — m^a^; 

Ci=migs+m2a3 + m^g2h + m^a2h-m^g2y—j ~ "^602 (^-^ J ; . (4.9) 



C2 = - fh\g2 - mla2] = -m\h + fhli-^\ . 



h> 

Here m*- denotes numerical parameters which have the form of some rational fractions. 
3. The time dependent elliptic 2-gap solutions, corresponding to the initial conditions 



3 and 4 in the system (|3.7] ) which are described by the expression (|4.l| ) at iV = 5, have 
the form 

5 

U{z) = 2^ p(z - ^,{t)) - 2(p(<^04) + P((^05)). (4.10) 

i=l 

In view of the initial conditions, the lower limits of the integration in ( [4. 3D have the common 
values </9oi|i<3 = and ip^^ = LOi, (/?05 = at the condition 3 and (/3o4 = —^^05 = ^2 at the 
condition 4. 

For computing the Aj-functions we substitute the expression ( 4.10| ) into the system 
( |3.5| ). Then, equating coefficients of the Laurent expansion in p of the left- and right-hand 
sides, we obtain the equation 

X^ + csX^ + C2X^ + ciX + Co = 0, 
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solutions of which describe functions n = 1,5. Coefficient functions Cn are de- 

scribed by the expressions 

5 

Co = ^(m°c/3 + n°a3)Fi(v9) + mlg2g3 + m^g?,a2 - ^§5203 

i=l 

5 

Cl = ^(m-g2 + n\a2)Fi{ip) - mlgj - m\g2a2 + m\al + mlaA\ 
1=1 

5 

C2 = mjgs + mlas; C3 = ^ mfFi{ip) - mlg2 + m|a2, 

i=l 

where 

Here mj, denotes some numerical rational fractions, f3i = 12p^((^) — g2, h = p{(p), 
ai = ai{p{^p)) denote some complicate functions, the explicit form of which we do not 
present here. 



4.3 Elliptic 3-gap solutions 

The possible types of the time dependent elliptic 3-gap solutions of the KdV equation 
(4.1) are determined by the values which are obtained from the compatibility condition 
between ( [l.lD and ( p. 9]) as t ^ 0. Under this condition, the number takes values 6, 8 in 
the formula (f4.1| ). 

1. The values = 6 and N = 7 determine two elliptic 3-gap solutions of the KdV 
equations with initial conditions 1 and 2 of the system (|3.9| ), which have the form 



Uiz,t) = 2Y,Piz-iPiit)), and t) = 2 ^ p(z - - 2p((^07), (4.11) 

1=1 1=1 

respectively. The time evolution of the poles (pi{t) are described by relations ( [4. 3] ) with the 
lower integration limits fo,i\i=YE = s-t the initial condition 1 and <Po,i\i=TE ~ ^' '^0,7 = 
at the initial condition 2. 

By substituting the expressions ( |4.11| ) into the finite-band equation ( |3.8D and equating 
coefficients of the Laurent expansion in p of the right-hand sides to zero, lead to the two 
equations 

6 7 

C6,»X* = 0, and J2 = 0, (4.12) 

i=0 i=0 



corresponding to the two 3-gap solutions with the values = 6 and N = 7, respectively. 
Coefficient functions cg^i and cy^j in ( 4.12 ) are rational functions on p{ipi) and p'{ipi). 
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Therefore Xi as roots of ( [4.12| ) are functions of ipi, i.e. Xi = Xi{ipi) where i = 1,6 and 
i = 1,7 at initial conditions 1 and 2, respectively. 

2. The value = 8 determines two similar elliptic 3- gap solutions of the KdV equations 
with the initial conditions 3 and 4 of the system (|3.9| ) , which have the form 



8 

U{z, t) = 2 ^ p(z - ^i{t)) - 2(p(^o7) + pM). (4.13) 

i=l 

The poles ^i{t) of ( [4.13| ) are described by relations ([4.3|), with the lower integration limits 
•^ai\i=Tfi = 0, (po7 = oji, ipo8 = ujj and (poili^j^ = 0, ipo7 = -f08 = -^3 at initial 
conditions 3 an 4, respectively. 

Substituting the expressions ( 4.13| ) into the finite-gap equation ( |3.8| ) and equating the 
coefficients of the Laurent expansion in p of the right-hand sides to zero, we obtain two 
equations of the form 

8 

Y,C6,^X' = 0, (4.14) 

i=0 

where Cj = Ci{ip). Eight solutions of (4.14) coincide with eight functions Xi{{pi), z = 1,8. 

The proposed approach is applicable for computing arbitrary elliptic n-gap solutions 
of the KdV equation. It can also be applied for computing finite-gap elliptic solutions for 
other integrable equations. 



5 Conclusion 

The solution of the KdV equation in the class of elliptic finite-gap functions is reduced 
to the solution of the system of finite-gap equations and auxiliary dynamic equations. In 
terms of the elliptic p-function this system reduces to simple algebraic relations which 
determin the parameters of the unknown solutions which are represented as rational func- 
tions of the p-function. This approach gives a simple algorithm for calculating arbitrary 
elliptic finite-gap solutions of the KdV equations at an initial time, which was demon- 
strated by the example of 1-, 2- and 3-gap solutions. The time evolutions of the unknown 
solutions with a known p-functional representation, are determined by the dynamics of 
their poles, which is described by coupled systems of dynamic equations. The latter always 
can be integrated with the help of the finite-gap equations. This was demonstrated by the 
example of elliptic 1-, 2- and 3-gap solutions. 

The above approach will also be applied to other integrable nonlinear equations in a 
future paper. 
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